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Abstract 

The CACTUS experiment recently observed a gamma ray excess above 50 GeV from the direction 
of the Draco dwarf spheroidal galaxy. Considering that Draco is dark matter dominated the gamma 
rays may be generated through dark matter annihilation in the Draco halo. In the framework of 
the minimal supersymmetric extension of the standard model we explore the parameter space to 
account for the gamma ray signals at CACTUS. We find that the neutralino mass is constrained 
to be approximately in the range between 100 GeV ~ 400 GeV and a sharp central cuspy of the 
dark halo profile in Draco is necessary to explain the CACTUS results. We then discuss further 
constraints on the supersymmetric parameter space by observations at the ground based ARGO 
detector. It is found that the parameter space can be strongly constrained by ARGO if no excess 
from Draco is observed above 100 GeV . 
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I. INTRODUCTION 



The existence of cosmological dark matter has been established by various astronomical 
observations. However, the evidences come mainly from the gravitational effects of the dark 
matter component. The nature of dark matter remains elusive and keeps one of the most 
outstanding puzzles in particle physics and cosmology Q]. The primordial nucleosynthesis 
and cosmic microwave background measurements constrain the baryon component and most 
dark matter component should be non-baryonic. The development in understanding the 
large scale structure formation requires the dark matter be cold. From the theoretical 
considerations the favored candidate for cold dark matter (CDM) seems to be the weakly 



interacting massive particles (WIMPs) |1|]. 

The WIMPs can be detected indirectly by observing the annihilation products, such 
as gamma rays, neutrinos, anti-protons and positrons. Exploring the anomalous results 
from the cosmic ray experiments is one viable way to identify the dark matter. Since the 
annihilation rate is proportional to the square of the dark matter density, the ideal sites 
for dark matter detection should have high dark matter density. The galactic center is 
believed to be a promising source of dark matter annihilation |2j. However, the existence 
of the central supermassive black hole and the supernova remnant Sgr A* contaminates the 
dark matter signals heavily . Alternative sites, such as the substructures of the Milky Way 
or the dark matter dominated dwarf spheroidal galaxies (dSph), have been studied in the 

RefsQ,y,y,y,Q. 

Recently, the CACTUS gamma-ray^ experiment reported an excess of gamma rays from 

Since Draco is dark matter dominated and no other 
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the excess has been attributed to the 
. The results are still preliminary and, 



the direction of Draco, a nearby dSph 
gamma ray sources are expected to be hosted 
annihilation of dark matter in the Draco halo 
if confirmed, will have important implications on the nature of dark matter and the density 
profile of Draco. Additional observations of the signal by other experiments is therefore very 
important. The GLAST a satellite-based experiment, and the MAGIC jl^J, a ground- 
based Atmospheric Cerenkov Telescope (ACT) , have been considered to check the CACTUS 
results [ulQ]. In the present work, we will discuss a possibility of detecting or constraining 
the gamma rays observed by CACTUS at ARGO [15j . a ground-based extensive air shower 
(EAS) detector. 
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In the next section we will first give the general formula for dark matter annihilation. 
Then we will discuss the implications of CACTUS results on the gamma ray spectrum and 
fluxes in Sec III. The sensitivity of ARGO is given in Sec. IV and the numerical results are 
presented in Sec. V. We conclude in Sec. VI. 



II. GAMMA RAYS FROM DARK MATTER ANNIHILATION 

The annihilation of two WIMPs can produce the continuous spectrum of gamma rays 
arising mainly in the decays of the neutral pions produced in the fragmentation processes 
initiated by the tree level final states. The fragmentation and decay processes can be simu- 
lated with the Pythia package [3]. 

The annihilation rate in unit time and unit volume is given by 

where a and v are the annihilation cross section and the relative velocity of the two dark 
matter particles respectively, n and p are the number and mass densities of dark matter and 
m is its mass, the factor 2 in the denominator arises due to the identical initial particles. 
We note that the annihilation rate is proportional to the square of the dark matter density 
and therefore, a high density region can greatly enhance the annihilation fluxes. 
The gamma ray flux from the Draco halo is therefore given by 

iy ! v y ' 2m? A-kD 2 4tt 2m 2 D 2 J An J H K 1 y ' 

where the halo profile is assumed approximately spherically symmetric with the density 
profile p(r), D = 75.8 ± 0.7 ± 5.4 kpc is the distance to Draco |l7(, 4> y (E) is the differential 
flux at energy E in a single annihilation in unit of 1 gamma GeV -1 . AQ represents the 
angular resolution of the detector. 

The density profile p(r) of Draco is constrained by observations. A recent analysis shows 
that both a cored and a cuspy profile, such as the NEW profi.e Q, are content with the 
observational data and the results of N-body simulation [19]. The 'astrophysical factor' in 
Eq. © defined as 

$ astro = Jp j AQ dn j 47rr WM , (3) 
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which is determined by the astrophysical quantities solely, is severely constrained by ob- 
servational data. It is found that <£> 'astro varies by a factor of approximately only 200, i.e., 
Castro = (3.2 x 10- 4 ~ 6.4 x lO" 2 ) GeV 2 cm" 6 kpc sr following Ref. Q- 

The other part in Eq. (J2J) is determined by particle physics which defines the nature of 
dark matter. We will calculate the 'particle factor' in the framework of the minimal super- 
symmetric standard model (MSSM). The MSSM is the most attractive model beyond the 
standard model of particle physics. In the R-parity conserved MSSM, the lightest supersym- 
metric particle, the lightest neutralino, provides a natural candidate for WIMP. The MSSM 
is well defined by a set of free parameters, which lead to the uncertainties in predicting 
the gamma ray flux from the particle physics. Once the particle factor is determined and 
combine with the astrophysical factor given above, we can give the predicted gamma ray 
flux from Draco. 



III. THE CACTUS EXPERIMENT 

CACTUS is a ground based Air Cherenkov Telescope (ACT) located at Solar Two near 
Barstow, California. CACTUS utilizes a set of 144 heliostats, each 42 m 2 , to form a com- 
posite mirror with a total effective area of about 6,000 m 2 . The threshold energy for gamma 
rays at CACTUS is about 50 GeV and the effective area for > 200GeV gamma rays reaches 
about 50,000 m 2 . 

Within the angular region of about 1° centered around the direction of Draco, CACTUS 
has recently observed an excess of approximately 30,000 photons for 7 hours observation 

n 

above the average background outside Draco |8|. The threshold energy of the photons is 
about 50 GeV. There is no significant excess observed if the cutoff energy is improved 
to about 150GeV. Although the results are still preliminary, however, if confirmed, the 
implications on dark matter are significant. It is interesting to consider the implications of 
the CACTUS experimental results seriously due to our completely ignorance of the nature 
of dark matter. In this section we will study the implications for the gamma ray spectrum 
and flux from the CACTUS results. 
The gamma events are given by 

r m x 

^observed = ^ / A eff (E)<$>(E)dEdQdT , (4) 
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FIG. 1: The spectrum of gamma rays from neutralino annihilation, with x = E^/m x , for the 
final state of W + W~, bb and rf. m x = lOO^OOGeV has been taken which gives almost identical 
spectrum for each final state. 

where e^n — 0.68 is the fraction of signal events within the angular resolution of the in- 
strument and the integration is for the energies above the threshold energy E t h and below 
the mass of neutralino, m x , within the angular resolution of the instrument AQ and for the 
observational time. The effective area A:// is a function of energy and <f>(E) = 0o^" is 
the flux of 7-rays from DM annihilation with (fi the intensity normalization and the 
shape of the spectrum. The effective area of CACTUS, which is energy dependent, can be 
parametrized as 

A eff ~ 47,000m 2 [i_e-°- 014 (^- 39 - 6GeV )] +11.9x E 7 (GeV) , (5) 

due to the simulation results [sj]. 

From Eq. (0J), we can see that in order to obtain the gamma ray flux from the observed 
event number we have to assume the gamma ray spectrum first. The spectrum of the gamma 
rays through neutralino annihilation depends on the final states into which the neutralinos 
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have annihilated. In Fig. ^we show the spectrum of gamma rays for the final states of gauge 
bosons XX ~ > W + W~ and for the final states of XX ~ > bb and XX ~^ r ^ which represent the 
two extreme cases that the annihilated gamma rays have soft and hard spectra respectively. 
In the figure we have plotted the spectrum for m x = 100, 500GeV respectively. We find the 
spectrum, expressed as a function of a dimensionless quantity x = E 1 /m x , is not sensitive 
to the mass of the neutralino, m x . 

The integrated gamma ray flux above the threshold energy of 50GeV is given by 

P m x r m x fjAf jy observed 

J 7 (> 50GeV) = / *{E)dE = / — fdE • / , (6) 

where we have assumed that the effective area has no zenith angle dependence within the 
angular resolution. From this equation we know the softer the spectrum the greater the 
gamma ray flux since A e ff is small at low energies. For soft spectrum, taking m x = lOOGeV 
and bb final states, we get Iy(> 50GeV) = 1.7 x 10~ 8 cm~ 2 s~ 1 , while for the hard spectrum, 
taking m x = 300GeV and rf final states, we get I 7 (> 50GeV) = 7.3 x 10 _9 cm _2 s~ 1 . This 
spectrum is taken in order not to give too much excess above 150 GeV. Concerning the 
uncertainties from the noise rejection procedures, the misidentification of the electronic and 
hadronic primary events and that the angular region of CACTUS is larger than that of Draco, 
the observed excess may be greatly larger than the real signal of dark matter annihilation. 
Therefore in our theoretical calculation we make an assumption that the uncertainty of the 
gamma ray flux is larger than the current CACTUS data by relaxing the lower bound by an 
order of magnitude. We finally get the gamma ray flux from Draco which is approximately 
in the range of 

7.3 x 10~ 10 < I 7 (> 50GeV) < 1.7 x lO^cm -2 *- 1 . (7) 

Since there is no significant excess observed above 150 GeV the gamma ray spectrum is 
further constrained. We assume that the events above 150 GeV do not exceed the Possion 
fluctuation of the background, which includes the misidentification of hadronic cosmic rays 
as gamma signals, the electronic comic ray events and the galactic diffuse gamma rays. We 
have adopted the expressions as 

4> h {E) = 1.49£7" 2 - 74 cm" 2 s" 1 sr~ 1 GeV'" 1 (8) 

for the hadronic contribution 

(j) e {E) = 6.9 x 10~ 2 E~ 3 - 3 cm- 2 s^sr- 1 GeV- 1 (9) 



for the electronic contribution 2l| . 



g aiac- 7 (^) = 8.56 x 10 _b E- 2 - 7 cm" 2 s _1 sr" 1 GeV'- 1 (10) 
for the Galactic 7-ray emission at the direction of Draco (I = 86.4°, b = 34.7°), extrapolated 



from the EGRET data at low energies 
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In principle the gamma ray flux above 150 GeV also depends on the spectrum of the 
gamma ray. However, due to Eq. (jHJ) the effective area above 150 GeV is not so sensitive 
to energy differing from that at energies below 100 GeV. Considering the large systematic 
uncertainties and the possible problems in the noise reduction procedure we approximate 
the effective area above 150 GeV as 50,000 m 2 , being a constant. Then we get a conservative 
upper limit of I 7 (> 150GeV). Assuming that about 90 percent of the hadronic comic ray 
background can be rejected within the angular region due to the different shape of the 
Cherenkov wavefront induced by electronic and hadronic showers, we get that 

I 7 (> 150GeV) < 3. x IO^otTV 1 . (11) 

In the next sections we will explore the supersymmetric (SUSY) parameter space to 
account for the gamma excess observed at CACTUS taking into account the constraints 
given by the Eqs. (J7|) and (JTTJ). 



IV. SENSITIVITY OF ARGO 



The ARGO-YBJ experiment, locates at YangBaJing (90.522° east, 30.102° north, 4300m 
a.s.l) in Tibet, China, is a ground-based telescope optimized for the detection of small size 
air showers. The energy threshold of the detector is designed to be about lOOGeV. The 
detector consists of a single layer of RPCs floored in a carpet structure covering an area of 
~ 10 4 m 2 . The detector is under construction and the central carpet will be completed early 
in 2006 and put in stable data taking soon after. 

The performances of the detector have been studied by means of Monte Carlo simulations 
. Defined as a product of the sampling area and the trigger efficiency, the effective area 
characterizes the power of the detector in recording the number of events for a given energy 
and time interval from a given direction. For both primary 7 and hadrons with energy near 
the threshold, the effective area can be approximately parameterized as A c s rs A W0 E 2A , 
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when the trigger condition is set to be larger than or equal to 20 fired pads, where Aiqq ~ 

100m 2 is the effective area for primary 7 ray events at the threshold energy of about 100 
I I 

GeV |23| . Above the threshold energy the effective area increases rapidly and reaches about 
10,000 m 2 for TeV gamma rays. At the same time, simulation also shows that at low energies 
the protons have lower trigger efficiency than photons. The effective area for protons near 
the threshold energy is about one order of magnitude smaller than that of gamma, leading 
to a great suppression of the background. 

The Draco dSph is within the field of view of the ARGO detector with the closest zenith 
angle to be ~ 27°. Following up observations on the gamma excess seen by CACTUS have 



been considered at GLAST and MAGIC 
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. Ground-based extensive air shower (EAS) 

n 

arrays with low energy threshold, such as ARGO |15| and the next generation all-sky high 
energy gamma-ray telescope HAWC j^J, have complementary properties to the satellite 
borne experiments and the ACTs. They have large effective areas and at the same time 
possess the advantages in large field of view and near 100% duty cycle. However, the EAS 
arrays usually have a poorer hadron-photon identification power. In this work, we will 
discuss how to constrain the gamma ray signal from Draco by the ARGO experiment. 

For this purpose, we focus on the events for the energy below ~ 400GeV, since we will 
see in the next section that the CACTUS excess constrains the neutralino mass to be lower 
than ~ AOOGeV. The number of background events for one year's data taking at ARGO is 
therefore also estimated in this energy range. To constrain the signal at the la level for one 
year's observation, the flux above 50 GeV from Draco is then constrained as 

IJ> 50GeV) = 2 4^- ^ , (12) 

^100^ e A cJ 100 (too j -jedh 

where again the zenith angle dependence of the effective area of the ARGO detector is 
ignored. 



V. NUMERICAL RESULTS 



In this section we will explore the SUSY parameter space to account for the CACTUS 
excess assuming that the excess (or a fraction of the excess) is generated by neutralino 
annihilation in the Draco halo. The constraint on the parameter space from ARGO is taken 
into account. 
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FIG. 2: The integrated 7-ray fluxes by neutralino annihilation from Draco above the threshold 
energy of 50 GeV as a function of the neutralino mass. The fluxes are given within the angular 
resolution of = 10~ 3 . Each point in the figure represents a set of low energy SUSY parameters 
which survive all the current limits. A boost factor 10 relative to the maximal astrophysical factor 
derived from [19| has been assumed. The lines shows the 2a constraints from the ARGO experiment 
assuming a W + W~ , bb or rf final state with or without gamma/hadron discrimination. 

The R-parity conserved MSSM is described by more than one hundred parameters de- 
scribing the soft supersymmetry breaking. However, for the processes related with dark 
matter production and annihilation, only several parameters are relevant under some sim- 
plifying assumptions, i.e., the higgsino mass parameter /1, the bino mass parameter Mi, the 
wino mass parameter M 2 , the mass of the CP-odd Higgs boson tua, the ratio of the Higgs 
vacuum expectation values tan/3, the scalar fermion mass parameter mj, the trilinear soft 
breaking parameter A t and To determine the low energy spectrum of the SUSY particles 
and coupling constants, the following assumptions have been made: all the sfermions have 
common soft-breaking mass parameters mr, all trilinear parameters are zero except those of 
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FIG. 3: Same as Fig. [21 except that a boost factor of 100 has been assumed. 

the third family; the gluino and wino have the mass relation, M3 = (a s (Mz) / o> em ) sin 2 Q\yM2, 
coming from the unification of the gaugino mass at the grand unification scale. However, 
to explore more general low energy phenomenological SUSY parameter space we relax the 
relationship between Mi and M2 derived from the grand unification scale. 

We perform a numerical random scanning in the 8-dimensional super symmetric param- 
eter space using the package DarkSUSY 26]. The ranges of the parameters are as follow- 
ing: 50GeV < \fi\, M 2 , M A , m f ~ < 5TeV, 1.1 < tan/5 < 60, -3m 9 - < A t , A b < 3m g ~, 
sign(/i) = ±1. The parameter space is constrained by the theoretical consistency require- 
ments, such as the correct symmetry breaking pattern, the neutralino being the LSP and so 
on. The accelerator data constrains the parameter further from the spectrum requirement, 
the invisible Z-boson width, the branching ratio of b — > S7 adopting the experimental data 
given by the Particle Data Group in the year of 2002 |27[. 

Another important constraint comes from cosmology. Combining the recent observation 
data on cosmic microwave background, large scale structure, supernova and data from HST 
Key Project the cosmological parameters are determined quite precisely. Especially, the 
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FIG. 4: Same as Fig. [21 except that a boost factor of 1000 has been assumed. 



abundance of the cold dark matter is given by |28( ^cdm^ 2 = 0.113' 



h0.008 
-0.009- 



We constrain 



the SUSY parameter space by requiring the relic abundance of neutralino < Q x h 2 < 
0.137, where the upper limit corresponds to the 3a upper bound from the cosmological 
observations. When the relic abundance of neutralino is smaller than a minimal value 
the thermally produced neutralino represents a subdominant dark matter component. We 
assume non-thermal mechanism to give the correct dark matter relic density [29]. The effect 
of coannihilation between the fermions is taken into account when calculating the relic 
density numerically. 

We find that a 'boost factor' at the order of 10 ~ 1000 is necessary to account for the 
CACTUS results. The 'boost factor' means that the astrophysical factor calculated by a 
cored or a cuspy profile in the Sec II should be enhanced by this factor to give the observed 
flux. The 'boost factor' requires a much sharper density profile compared with the NFW 
profile, such as a Moore profile [30 ] or a spike profile due to the existence of an intermediate 
mass central black hole |3l( in Draco. 

In Figs. |2l El and |H we plot the integrated 7-ray fluxes above the threshold energy 50GeV 
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within the solid angle AQ = 10~ 3 as a function of the neutralino mass. The results in Figs. 
El El and El have enhanced the astrophysical factor by a boost factor of 10, 100 and 1000 
respectively. Each point in the figure corresponds to a model with a set of definite SUSY 
parameters in the 8-dimensional parameter space which can explain the CACTUS results 
constrained by Eqs. (J7J) and (|11|) and allowed by all other collider and cosmology constraints. 
The scatter of the points represents the uncertainty coming from the unknown soft SUSY 
breaking parameters. The lines represent the 2a constraints of ARGO by assuming the 
final states being the W + W~ , bb or rf. For the upper set of line we have assumed no 
hadron/photon discrimination at all, while for the lower set of lines we assume a part of 
hadrons are rejected based on neural network so that the significance of detection is improved 
by a quality factor of 1.6 3^|. From the figure we can see that if no excess is observed at 
ARGO above 100 GeV a large part of the parameter space is constrained. 

It is worthwhile commenting on the results here. First, if extending the gamma spectrum 
to lower energies, we find the CACTUS result is difficult to reconcile with the EGRET result 
j^J which did not observe excess at the direction of Draco between 1~ 10 GeV. Therefore 
a hard spectrum is expected to reconcile the EGRET and the CACTUS results, which 
requires the dominant annihilation product be rf ll|. The hard spectrum leads to more 
opportunities to observe the signal in ARGO which can be seen from the Figs. El and 0] 
Alternatively one would assume that only about 1 percent of the present excess is real signal 
from the annihilation of the dark matter. In this case we find parameter space to account for 
the signal and be consistent with EGRET result in the range of 250GeV < m x < 800GeV. 
The parameter space can be constrained by ARGO only for the rf final states. Second, the 
CACTUS result may also imply a monochromatic gamma spectrum at the energy of about 
50 GeV. However, it is found that the branching ratio that two neutralino annihilate into 
two photons should be more than a half to be consistent with the EGRET result, which 
is incompatible with the SUSY model 3|. Finally, if we assume that only about ^ 1% 
of the excess comes from DM annihilation, the signal can be explained without introduce 
any 'boost factor' if taking the nonthermal mechanism into account. This may be a natural 
assumption, while the confirmation of the gamma events from DM annihilation requires an 
instrument with better angular resolution, such as GLAST [3| to suppress the background. 
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VI. SUMMARY AND CONCLUSION 



In this paper we have discussed the possibility of constraining the signal observed by the 
CACTUS experiment at the ground based EAS detector, ARGO. We assume the excess of 
gamma rays observed at CACTUS is produced by supersymmetric dark matter annihilation. 
We then explore the SUSY parameter space to give signal consistent with the CACTUS result 
and discuss the possibility to constrain the parameter space at ARGO. Our calculation shows 
that, depending on the gamma spectrum, ARGO will be able to constrain a large part of 
the parameter space if no signal is detected for one year observation. 

If the CACTUS signal is finally confirmed, the implication on dark matter is dramatic. 
The central cusp of the dark halo at Draco should be much sharper than that of a NFW 
profile. The neutralino mass should be at the range of 100 ~ 400 GeV to explain the signal 
of CACTUS. Furthermore, the spectrum of the annihilation gamma ray should be very hard 
in order to be consistent with the EGRET null result at the direction of Draco at the energy 
range between 1 GeV and 10 GeV. 
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